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Abstract. In this paper, I have attempted to highlight key results from deep 
extragalactic surveys at mid-infrared wavelengths. I discuss advances in our 
O , understanding of dust enshrouded star-formation and AGN activity at < z < 3 

from IRAS, ISO and Spitzer. The data seem to indicate that about 70% of 
the co-moving star-formation rate density at 0.5< z < 3 is obscured by dust 
and that AGN, including obscured sources, account for <20% of the co-moving 
■ bolometric luminosity density. There is tentative evidence that the mode of 

star-formation changes as a function of redshift; star-formation at z ~ 2 is 
preferentially in massive, ultraluminous infrared galaxies (ULIRGs) while z ~ 1 
^ ' sources are luminous infrared galaxies (LIRGs) which are about 1 mag fainter 

| than ULIRGs in the near-infrared. This evolution of the star-formation mode, is 

CO . similar to the evolution of the redshift distribution of X-ray sources as a function 

of X-ray luminosity and would suggest an extension of the downsizing hypothesis 
C^l | to include both AGN and star-forming galaxies. Measuring dust-enshrouded 

star-formation at z > 3 will become possible only with future facilities like 
ALMA. Currently the presence of dust can only be assessed in a small fraction 
of the youngest starbursts at z > 5 by looking for redshifted large equivalent 
f~| . width Ha emission in broadband filters like the IRAC 4.5/im passband. Ha 

Oh" to UV ratios in these objects are a tracer of dust extinction and measuring 

1. \ this ratio in GOODS galaxies indicate dust in ^20% of star-forming galaxies at 

z > 5. Finally, implications for reionization based on the measured stellar mass 
' density and star-formation rates of galaxies at these redshifts are discussed. 
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1. Pros and Cons of the Mid-infrared 

The mid-infrared regime of the electromagnetic spectrum, is generally defined by 
the wavelength range which can be detected with arsenic/antimony doped silicon 
semiconductors and spans 5— 40/xm. Extragalactic surveys in the mid-infrared 
are typically limited by background noise since the ther mal emission from zo - 
diacal light in the solar system peaks between 10— 20/um (jKelsall et al. 1998). 



These wavelengths correspond to the trough between the optical/near-infrared 
and far-infrared wavelengths in the spectrum of the extragalactic background 
light (EBL). To illustrate, the EBL intensity derived by integrating the light 
from galaxies detected at 15 and 24//m is around 2.5-5 nW m -2 sr _1 while the 
EBL intensity at both 140 //m and 2.2 fim has a value of about ~25 nW mT 2 sx~ l 
(|Lagache et al. II2005I L Since the EBL is the sum total of redshifted emission 



from sources integrated over the observable Universe, the bulk of which arises 
from galaxies at z < 1.5, this implies that the energy emitted at near- and 
far-infrared wavelengths is almost an order of magnitude larger than that at 
mid-infrared wavelengths. This is generally true even for individual star- forming 
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galaxies; the mid- infrared regime lies between the optical/near-infrared stellar 
photospheric emission and the far-infrared emission from dust heated to T~20- 
40K, both of which together dominate the bolometric luminosity of the galaxy. 

The mid-infrared is a complicated spectral region representing the inte- 
grated emission from broad poly cyclic aromatic hydrocarbon features, contin- 
uum emission from dust grains between 300K and 40K and silicate absorption 
dDesert et al. 1119901 ). As a result, even small errors in redshifts translate into 
large errors in spectral k— corrections making estimates of the true rest-frame 
luminosity difficult (Figure 1). 

Despite these difficulties, the mid-infrared is the most powerful tracer of 
dust enshrouded star-formation activity out to redshifts of 3 (Figure 2). This is 
possible because of the tight correlation b etween mid- and far- i nfrared emission 
on b oth sub- galactic and galactic scales ( Calzetti et al. 20051 ; ICharv &: Elbaz 



l200ll ). For example, rest-frame 7/xm emission has been shown to trace Ltr with 
a la scatter of 40% in the local Universe. The correlation appears to hold 
even in the z ~ 1 Universe as determined from a comparison between 7fim, 
70fim and 1.4 GHz radio luminosities of individual galaxies. However, there 
is some evidence that at lo w metallicities (Z<l /4 Zq) the PAH line strengths 
show significant differences ( Smith et al. 20061 ) . On the instrumentation side, 



mid-infrared observations have arcsec sized beams even with small ~lm class 
telescopes. This has the benefit of alleviating the effects of source confusion. It 
also simplifies the identification of the counterparts of mid-infrared sources at 
optical/near- infrared wavelengths and thereby enables the measurement of their 
spectroscopic redshifts. Finally, mid-infrared arrays suffer less from latents and 
transient effects in their responsivity than far-infrared arrays and need to be 
cooled only down to a few degrees Kelvin instead of the milli-kelvin requirements 
of submillimeter bolometers. 



2. Key Extragalactic Mid-infrared Surveys 

Deep IRAS 12 fim and 60 //m observations in the direction of the North Ecliptic 
Pole ([Hacking fc Houck II1987T ). were the first to provide evidence for redshift 



evolution of the infrared luminosity function of galaxies. The survey measured a 
source density of ~30 deg -2 down to 10 mJy which was a factor of 1000 deeper 
than the IRAS All Sky Survey. The source counts down to this flux density limit 
provided tentati ye evidence for evolu tion of the infrared luminosity function to 
redshifts of 0.2 ([Hacking et al. Ill987ft . 



Subsequent mid-infrared surveys which have revolutionized our understand- 
ing of galaxy evolution at c osmological distances are those undertaken with 
ISOCAM at 15 /mi surveys (|Genzel & Cesarskv I l2000h and Spitzer /MIPS at 



24/im. 

ISOCAM undertook a variety of deep 15/im surveys, the deepest of which, 
apart from those in lensing cluster fields, was in the Hubble Deep Field North 
( Aussel et al. 19991 ). These observations which covered an area of ~ 20 arcmin 2 



centered on the HDF-N, were sensitive down to a flux limit of 100//Jy and re- 
liably detected about ~40 sources. The differential source counts was more 
than an order of magnitude larger than expected if the local infrared luminosity 
function of galaxies did not evolve with redshift. This data provided reliable 
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evidence for a strong ~(l+z) 4 luminosity evolution of infrared luminous galax- 
ies. Identification of the optical counterparts of the 15/iin sources and measure- 



ment of their redsh i fts fr om ground-based spectroscopic surveys (jWirth et al 



20041 : ICohen et al~l |2000| ) revealed that these were galaxies at z ~ 0.5 — 1.0. 



Cross-correlation with X-ray surveys from Chandra and XMM (Brandt, these 
proceedings) as well as their optical spectra revealed that only about 20% of 
these objects appear to have AGN in them implying that star- formation was 
powering the bulk of the emission. As a result, one could apply the local cor- 
relation between mid- and far- infrared luminosities, which revealed that these 
were infrared luminous galaxies with L/r=L(8— 1000/xm)> 10 11 Lq with the 
majority of them being LIRGs (10 11 < L/^j < 10 12 Lq). The estimates of their 
star-formation rates are in the range of 10-100 Mq yr _1 , most of which is dust 
obscured. For the subset of sources which are radio detected, the derived far- 
infrared luminosities of these sources agree well with the radio-FIR correlation 
in the local Universe providing evidence against an evolution of the far-infrared 
spectral energy distribution with redshift. The SED fitting also yielded a contri- 
bution of ISOCAM galaxies to the DIRBE observed far-infrared background of 
65±30% (El baz et al. 20021 ). The contribution of dusty starbursts to the total 



EBL has since been confirmed more recently by s tacking Spitzer far-infrared 
data on mid-infrared sources bv lDole et al~l (|2006l ^. 

Phenomenological models which attempted to fit the source counts from sur- 
veys at various wavelengths, based on local galaxy templates, i ndicated a strong 



evolu t ion of the far-infrared lu minosity density with redshift (jCharv &: Elbaz 



l200ll ; iFranceschini et al. 200 ll ). These models revealed that the far- infrared 
luminosity density is between a factor of 5—10 higher than the UV luminosity 
density at z ~ 1 implying an increase in the average dust extinction in the 
Universe from E(B-V)~0.15 at z ~ 1 to E(B-V)~0.3 at z ~ 1. However, it is 
difficult to have these models account for cosmic variance which cause source 
counts to vary by as much as a factor of 2 between different fields. Also, phe- 
nomenological models are degenerate in factoring in the amount of luminosity 
evolution and density evolution in the infrared luminosity function due to the 
difficulty in measuring the faint end slope of the luminosity function. Various 
models which simultaneously fit the source counts at different wavelengths have 
slightly different parameters for the redshift evolution of LIRGs and ULIRGs 
( Chary &: Elbaz Il200ll ; iLagache et al. 1 12004 ) . In addition, the underlying mid 
and far-infrared templates used in the models strongly affect the evolution pa- 
rameters. Nevertheless, the evolution of the far-infrared luminosity density ap- 
pears to be robust to these parameters. 

There were two significant problems. Due to the fact that the PAH features 
get redshifted out of the ISOCAM 15//m bandpass at z ~ 1.2, it was impossi- 
ble to assess evidence for evolution of the infrared luminosity function at higher 
redshift without using galaxy counts at 850/um which are discussed elsewhere 
(R. Ivison, these proceedings). Secondly, the integral of the co-moving star- 
formation rate density from the models overpredict ed the measured comoving 
stellar mass density seen bv lDickinson et al. ( 20031 ). This could either be due 



to the stellar mass density being underestimated from the uncertain contribu- 
tion of the high Mass/Light ratio old stellar component in galaxies. Alternately, 
it could be due to the fact that at high redshifts, the mid-infrared determined 
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bolometric luminosity densities are overestimates. Deep, wide Spitzer surveys 
(e.g. Dickinson et al. Spitzer GO-3) are well positioned to address this prob- 
lem by providing 70pm coverage of large numbers of z > 2 galaxies, enabling 
bolometric luminosities to be validated using stacking analysis. 



3. GOODS MIPS 24pm Survey 

The Great Observatories Origins Deep Survey 24pm observation s were taken 
with the Multiband Imaging Photometer Spectrometer (MIPS; iRieke et al. 



(|2004T n. The 24pm passband can detect PAH emission out to z ~ 3. How- 



ever, unless the observations are adequately deep, observations would be un- 
able to detect typical infrared luminous galaxies much beyond the limits set by 
ISOCAM. For example, at a flux limit of 80/xJy, 24 micron observations are 
sensitive to a typical 3xlO n Lq out to z ~ 1.2 which is about the same as what 
the ISOCAM HDF observations could detect such objects out to. By going 
down to 20/xJy, such galaxies are now being detected out to redshifts 2.5, en- 
abling a comparison with the z ~ 2 submillimeter galaxy population and z ~ 3 
Lyman-break galaxy population. 

The Spitzer beam size at 24pm is 5.7" FWHM which raises the issue of 
source confusion. However, due to the accurate pointing of Spitzer and ~0.2" 
rms with which we can align the 24pm images with observations at 3— 8pm, 
the concerns of confusion noise placing a floor on the sensitivity are misplaced. 
The surface density of sources is 17 arcmin -2 in GOODS-N and 24 arcmin~ 2 in 
GOODS-S. Also, at this resolution, the vast majority of extragalactic sources 
are point sources. Using a prior position based point source fitting algorithm, 
we have been able to extract sources with a completeness of 84% at 24/xJy 
(Figure 3). The 5a point source sensitivity limit of this survey is 25/iJy. The 
fundamental reason for the success of this technique is because every MIPS 
source is detected in deep 3.6^m observations. 

The differential 24^m source counts are similar in shape to the 15/im counts 
with a break in dN/dS at 400/uJy (Figure 3). The slope of the counts at fainter 
fluxes is ~S~ 16 . One of the most salient features of the 24pm counts is they peak 
(in dN/dSxS 2,5 ) at fainter flux values than predicted from the phenomenological 
models which are fit to the multiwavelength counts. Attempts to resolve this 
have included chan g ing th e evolution of the LIRG to ULIRG ratio with redshift 



as m 



-.Chary et al. I (12 004) or weakening the strength of the 11 and 12pm PAH 
(jLaeache et al. 1120041 )" in the galaxy templates so that objects which contribute 



to the peak of the Ihpm counts contribute at fainter fluxes to the 24pm counts. 

There is little observational evidence for a weakening of the 11 and 12^m 
PAH (e.g. Pope et al. 2006, these proceedings). Stacking of IRS spectra of 24pm 
sources indicate that apart from an underlying hot dust continuum which might 
lower the equivalent widths of all PAH between 6 and 12pm, the relative line 
ratios do not vary except for metal poor galaxies. Since dusty starbursts appear 
to be in massive galaxies, which have also enriched the ISM in the process of 
increasing the dust content, there is less concern that the bolometric luminosities 
for the bulk of the galaxies are incorrect although some recalibration of the 
libraries probably have to be done for the SEDs of the most luminous galaxies 
(Ljr > 10 13 Lq) for which no low-redshift counterparts exist. 
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4. Evidence for "Downsizing" in AGN and Galaxies 



From the redshifts and 24/xm flux densities of individual sources, we can calculate 
their far-infrared luminosity using the template fitting approach ( Chary &: Elbazl 
12001] ) . The accuracy of this approach has been verified at z < 2 from 70/iin 
and 850/um stacking of LIRGs and ULIRGs. Also, the spectral energy distri- 
bution of galaxies at z ~ 2 selected using the BzK tech nique also appear t o 
match the templates spectral energy distribution closely ( Daddi et al. 1 I20051 ) . 
The far-infrared lumino sity can be converted in to a star-formation rate using 
the conversion factor in iKennicutt et al. ( 19981 ). One possible uncertainty in 



this conversion is estimating the contribution to the far-infrared luminosity from 
dust being heated by the evolved stellar population. Deep optical surveys pro- 
vide a measure of the star-formation seen in the rest-frame ultraviolet. Since the 
bulk of the far-infrared emission is dust heated by the ultraviolet radiation field 
which is dominated by star-formation, adding the UV and FIR star-formation 
rates yields the star- formation rates of individual objects. 

Although the stellar masses of galaxies can be derived using SED fits to the 
multiband photometry, a comparison between the 3.6/mi magnitudes of LIRGs 
and ULIRGs which traces rest-frame near-infrared emission provides a good 
estimate of the stellar mass. The photometry indicates that LIRGs at z ~ 1.5 
have a median 3.6^m brightness of 21.2 AB mag while ULIRGs are 20.5 AB 
mag. At z ~ 2.4, LIRGs are 22.5 AB mag while ULIRGs are 21.5 AB mag. 
Thus ULIRGs are systematically, about 1 AB mag brighter at 3.6/xm indicating 
that they are more massive. 

Finally, AGN can be estimated in the 0.5 < z < 3 range either by identi- 
fying power laws from hot dust in the observed 5.8— 24/xm range. Alternately, 
the X-ray luminosity and photon index of sources provides a good measure of 
AGN activity. We caution that Compton-thick AGN are undetected even in the 
deepest Chandra surveys beyond z ~ 0.6 and the IRAC power law identification 
of AGN breaks down at redshifts of z > 2 because dust sublimes at around 
~1500— 2000 K corresponding to rest-frame ~2/xm. At z ~ 2, the 2/im emitting 
hot dust is redshifted out of the 5.8/im bandpass. 

Putting these various datasets together, we confirm the ISOCAM result im- 



plying an evolution of the far- i nfrared luminosity density w ith redshift (jBabbedge et al 



20061 ; iLe Floc'h et al. Il2005l : iPerez-Gonzalez et al. I [2005). However, the data 



also indicate an evolution of the galaxies dominating the co-moving star-formation 
rate density. At z ~ 1, LIRGs which typically have stellar masses of ~10 9 ~ 10 Mq 
and star-formation rates of 10-100 Mq yr _1 are dominating the energetics. At 
z ~ 2, it is the ULIRGs which have stellar masses >10 10 Mq that are dom- 
inating the star- formation rate density (Figure 4). The AGN contribution to 
the luminosity density at both redshifts appear to less than ~20%. Thus, lower 
mass galaxies form later in cosmic time and more massive galaxies build their 
stellar mass earlier. 

This trend is similar to the redshift dist ribution of X-ray sources as a func- 
tion of X-ray luminosity ( Ueda et al. 1120031 ). The X-ray data seem to suggest 



that the most massive (~10 9 Mq) nuclear black holes built up their mass in a 
luminous accretion event around z ~ 2 — 3 while the ULIRGs were building up 
their stellar mass. The lower mass black holes on the other hand, grew around 
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z ~ 1 when the LIRGs were building up their stellar mass providing impetus 
for the downsizing scenario extending to AGN and starbursts. Since low mass 
galaxies are more numerous and populate the Hubble sequence, this implies that 
the Hubble sequence was probably in place only at z < 1. 

However, there are potential systematics. There is some evidence in our 
IRS spectra that the color temperatures of the mid-infrared spectra are rela- 
tively cool. Templates of local galaxies which are used to derive the bolometric 
corrections typically have steeply rising hot dust continuum with increasing L/^. 
This would imply that bolometric corrections might have been overestimated and 
many of the objects are less luminous. The effect is particularly important in 
the 2 < z < 3 regime where we have limited checks on the shape of the template 
from stacking analysis. There is also evidence in some objects that continuum 
from an AGN results in the 6.2/mi PAH equivalent width being reduced by a 
factor of 3 compared to local starburst galaxies such as M82. This would again 
imply an incorrect bolometric correction in the z > 2 sources. Finally, due to 
the difficulty in identifying AGN at z > 2 from the IRAC colors, we might 
be underestimating the AGN contribution which would be predominantly to 
the brightest 24^m sources thereby lowering the ULIRG contribution. Further 
checks, especially in the 2 < z < 3 range, where both photometric and spec- 
troscopic redshifts are limited are essential to verify the result shown in Figure 
4. 



5. Dust at Higher Redshifts and Implications for Reionization 

Detecting dust emission at redshifts greater than 3 is difficult. The strongest 
PAH features at 6.2 and 7.7 fiva. are redshifted out of the MIPS 24^m band- 
pass. Only source with large AGN contributions displaying significant hot dust 
emission are detected at z > 3 at 24//m. It appears that the equivalent width of 
the 3.3/mi PAH is low enough that it would be undetectable above the stellar 
continuum and the correlation between 3.3/xm emission and Ljr is not well es- 
tablished. 70//m surveys are about one to two orders of magnitude less sensitive 
than required to detect dust in all but the most hyperluminous sources. As a 
result, at z > 3, 850/im surveys are best suited for detecting dust emission but 
the current sensitivity limits (Figure 2) will not be significantly improved until 
the Atacama Large Millimeter Array (ALMA) becomes operational. 

One possible technique has been to measure redshifted Ha emission through 
broadband phot ometry. This techn ique can only be applied over a limited red- 
shift range. In ICharv et al.l ( 20051 ). we demonstrated the possibility that the 



enhanced 4.5/xm/3.6/im flux ratio of the lensed z = 6.56 galaxy HCM6A is 
due to large equivalent width Ha. A study of the broadband photometry in 
GOODS galaxies has since revealed that about 10—20% of spectroscopically 
co nfirmed z > 5 obje cts show a similar enhancement in the 4.5/mi flux (Figure 
5; Chary et al. 2X106). The Ha determined star- formation rates are on average 



higher by a factor of ~3, compared to the UV determined star-formation rate, 
indicating the presence of dust in a small fraction of high redshift objects. The 
presence of dust within 1 Gyr of the Big Bang argues for the production of dust 
from SNe since red giants and AGB stars have had insufficient time to evolve. 
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Since the dust content correlates with the number of SNe, it must increase 
with decreasing redshift and is presumably directly correlated to the increase 
in stellar mass. This implies that for a small fraction of objects, the measured 
star-formation rate from the UV continuum is lower than the true star-formation 
rate. Since typical starbursts have an e-folding time associated with them, the 
conclusion is that that the star- formation rate (p*(t)) of a galaxy at some time 
t < t' in its past history must be: 

P*(t) = p*(t') exp[-(t - t')/T SFR ] exp[-r(t) + r(t')] (1) 

where, 

r(i') oc / p*(t)dt (2) 



o 

r SFR is the e-folding time for the star-burst while r(t) is the extinction 
in the star-formation rate metric (i.e. UV wavelengths) at time t. Naturally, 
r(t)<r(f). 

We have now been able to measure the star-formation rate and stellar mass 
density in galaxies at z ~ 6. For canonical estimates of the escape fraction 
(fesc~0.15) and clumping factor (C=30), the measured star-formation rate at 
z > 6 falls more than a factor of 10 short of the minimum needed to reionize 
the IGM. However, if the dust content is increasing with redshift, the past star- 
formation rate must be higher not just by the e-folding time of the burst but also 
by the increase in dust extinction which is directly proportional to the stellar 
mass. 

The past history of star-formation can either be identified by fitting SEDs to 
the multiband photometry and extrapolating backwards in time or by fitting an 
evolving starburst model to the measured star-formation rate and stellar mass 
density at z ~ 6 (Figu re 6). From rest-fram e ultraviolet observations of Lyman- 
break galaxy samples ( Bouwens et al. 20061 ). it appears that the star- formation 



rate density is decreasing at z > 6, while from the extrapolation process the 
star-formation rate density must be increasing at higher redshifts to reproduce 
the observed stellar mass density. This discrepancy is probably attributable to 
the large completeness corrections associated with identifying galaxies at z > 7 
using current instrumentation. 

The primary conclusions are as follows: 1) If the measur ed star-formation 
rate density and s t ellar mass density at z ~ 6 are acc urate ( Yan et al. 1 1 20061 : 



Chary et al. 1120061 : lEvles et al.l 12005 : IStark et al.ll2006l ). then reionization must 



be a brief, inhomogeneous processing lasting only 20-50 Myr at 6 < z < 7; 
2) If the process of reionization starts at z ~ 9 and is prolonged, then the 
stars that dominate the reionization process must not contribute to the stellar 
mass density at z ~ 6 i.e. the mass of stars and remnants at z ~ 6 must be 
higher than the measured stellar mass by a factor of ~10. This would suggest 
a contribution from high mass, zero metallicity Population III stars or from 
the low mass end of the galaxy mass function. Either of these scenarios are 
unnecessary if the escape fraction and clumping factor are significantly different 
from the canonical estimates of 0.15 and 30 respectively. 
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Figure 1. Uncertainty in determining Lir from monochromatic observed 
24/im fluxes due to an error in the photometric rcdshift estimate. Due to the 
interplay of emission and absorption features in the mid-infrared, small errors 
in redshift estimates result in the incorrect far-infrared template being chosen 
from template libraries. This manifests itself into large errors in the bolomet- 
ric correction and thereby the derived bolometric luminosity of objects. The 
la scatter in the derived Lir due to typical optical/near- infrared photometric 
redshifts is a factor of 3 for a la scatter in Az of 0.24 which corresponds to a 
cr(Az/(l+z spec )~0.1. 




Redshift 



Figure 2. Sensitivity of deep surveys at different wavelengths to 
Lir=L(8— 1000 /im), a proxy for star-formation rate. Also plotted are sources 
with either photometric or spectroscopic redshifts. Stars indicate X-ray 
sources i.e. AGN, while circles are star-forming galaxies. The GOODS 24 /xm 
observations are at least an order of magnitude more sensitive to dust ob- 
scured star-formation and AGN activity than other surveys. 
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Figure 3. Differential cxtragalactic source counts at 24^m in the GOODS- 
N field. The empty squares are the measured counts, the solid squares are the 
counts corrected for completeness, the open diamonds are the completenes s 
corrected counts in the GOODS test field in ELAIS-N 1 (ICharv et al. 1120041) 
while the asterisks are the counts in various other fields (|Papovich et al. ! 2004; 
iMarleau et al. I r2004f ) . Application of the prior-based deblending technique 
results in an extracted source catalog which is 84% complete at 24^iJy illus- 
trating that source confusion is not a significant source of noise in GOODS. 
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Figure 4. Relative contribution of LIRGs (hatched region), ULIRGs (back- 
hatched region) and AGN (stars) to the co-moving far-infrared luminosity 
density between < z < 3. The dotted symbols underlying are various mea- 
sures of star- format ion from the UV and Ha referenced in ICharv fc Elbaz I 
(2001). The systematic uncertainty in this measurement is a factor of 3 — 4 at 
z > 2 due to the limited number of sources with spectroscopic redshifts and 
uncertainty in the AGN contribution (See text). 
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Figure 5. Identification of Ha in emission using bro adband colors. Cir- 
cled sources are c andidate Ha emitters with z spec > 5 (|Charv et al. Il2006t 
ICharv et al.l l2005) . Comparison of the Ha SFR to the UV continuum derived 
star-formation rate provides a unique technique to measure dust obscured 
star-formation at z > 5. However, this technique can be used only in the 
most violently star-forming galaxies because the equivalent width of Ha must 
be very high for it to be detectable in a broadband filter like Spitzer /IRAC 
4.5/im. 
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Figure 6. Figure 1: (Left Panel) The solid black line shows the minimum 
star-formation rate as a function of redshift, for a Salpeter IMF, which is 
required to reionize the IGM (escape fraction of 0.15 and clumping factor of 
30). The black symbol shows the measured star- formation rate density at z^6 
(Bouwens et al. 2006). The solid grey line, shows the star- formation history 
of the Universe required to reproduce the measured stellar mass density as 
well as the star- formation rate density at z^6. Since the grey line can exceed 
the solid line for a brief interval of time, reionization must be a brief process 
lasting 20-50 Myr at z~6.5. (Right Panel) If on the other hand reionization 
took place at z~9, then either the stellar mass density at redshift 6 must be 
xlO higher than what is measured due to a steep faint end mass function 
or the sources responsible for reionization at z^9 do not contribute to the 
stellar mass density at z^6. The dashed grey lines show the role of increasing 
extinction on the measured star-formation rate density for Ay=0.5 mag and 
Ay of 1 mag. 



